Intertwining of exactly solvable generalized Schrodinger equations 



(N 



A.A. Suzko"'^ E.P.Velicheva" 

Joint Institute for Nuclear Research, 141980 Dubna, Russia and 
^ JIPENR, National Academy of Sciences of Belarus, Minsk 

The Darboux transformation operator technique in differential and integral forms is applied to 
the generalized Schrodinger equation with a position-dependent effective mass and with linearly 
energy-dependent potentials. Intertwining operators are obtained in an explicit form and used 
for constructing generalized Darboux transformations of an arbitrary order. A relation between 
supersymmetry and the generalized Darboux transformation is considered. The method is applied 
to generation of isospectral potentials with additional or removal bound states or construction of 
new partner potentials without changing the spectrum, i.e. fully isospectral potentials. The method 
is illustrated by some examples. 
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I. INTRODUCTION 



The problem of exact solvability of the Schrodinger equation has been extensively considered since the beginning of 
^ quantum mechanics [l| . A factorization technique was introduced by Schrodinger and was used to solve the harmonic 
: oscillator, the hydrogen atorn, the Kepler motion. It should be noted that factorization is closely connected with 
I ' the Darboux transformation [2*1 and the supersymmetry approach introduced by Witten in quantum mechanics 
, Many interesting exactly solvable models have been constructed for Schrodinger equation using the Darboux and 
• Bargmann transformation techniques [1]- 14|. Many advances have been made in the field of classifying quantum 
^ ' mechanical potentials according to their symmetry properties and in the area of different applications of ex actly 
solvable models in quantum mechanics, in particular, in atomic and nuclear physics, in statistical mechanics (7j-[l4|. 
I Potentials of Bargmann- type and corresponding exact solutions can be obtained from the integral equations of the 
I r inverse scattering problem in the Gelfand-Levitan and Marchenko approaches [l^ with the degenerate kernels 
^ of the transformation operator. The differential Darboux transformations (or the method of itertwining) are closely 
, related to supersymmetry and have a variety of applications in different fields of physics. The technique of iterated 
' integral operators for Schrodinger equation is based on the integral Gelfand-Levitan and Marchenko equations (see 
^.1. '. e.g. The differential and integral approaches of obtaining exactly solvable equations are very close to each 

• ' other, but sometimes the differential method is more convenient and sometimes the integral method is preferable. In 
particular, the integral form of transformations is used for construction of phase-equivalent potentials. Suggested in 
L.J . [l2l [l3| generalized Darboux and Bargmann transformations permit construction of exactly solvable potentials for 
I ' variable values of energy £ and angular momentum I. The transformations in differential and integral forms were 
• • , considered and it was shown how to construct the phase-equivalent potentials. Later on, this approach was generalized 
. ^ on the Schrodinger equation with weighted energy [3] and recently [l^l the intertwining technique has been applied 
k> \ to construct exactly solvable potentials for this equation. The Schrodinge r eq uation with energy- dependent potentials 
5^ is widely used in nuclear physics, helium clusters and metal clusters [21|-[25|. 

■ In the last few years the research efforts on the topic of algebraic methods have been considerably intensified [26|-[36j 
due to the rapid development of nanoelectronics, the basic elements of which are low-dimensional structures such as 
quantum wells, wires, dots and superlattices (stI. Issj. For investigation of nonuniform semiconductors, in which the 
carrier effective mass depends on position, the generalized Schrodinger equation with position-dependent effective mass 
is used [39.] . One of the most important problems of quantum engineering is the construction of multi-quantum well 
structures possessing desirable spectral properties. The technique of intertwining relations or Darboux transformations 
allows one to model quantum well potentials with a given spectrum. Although Darboux transformations are applicable 
to the position-dependent mass Schrodinger equation [28, 2SL 34, iiS] , the Schrodinger equation with weighted energy 
[20j and the generalized Schrodinger equation with position-dependent mass and weighted energy [H, [sl], however, 
their applications are very complicated for realization and there remain many problems to study. Nowadays there 
is a great interest to generate Hamiltonians with a prescribed energy spectrum. However, these problems are worth 
investigated for all cases of the generalized Schrodinger equations except for the standard Schrodinger equation ( see 
e.g. [5]-[8]). 

In this paper we will focus on Darboux transformations of an arbitrary order applied to the generalized Schrodinger 
equation with a position-dependent mass and energy-dependent potentials. These generalized Darboux transforma- 
tions include all known cases: Darboux transformations for the Schrodinger equation with the position-dependent 
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mass, the case with hnearly energy- dependent potentials as weh as Darboux transformations for the conventional 
Schrodinger equation. It is known that Darboux transformations and supersymmetry are recovered for the standard 
Schrodinger equation. Here, we consider interrelations between supersymmetry and Darboux transformations for 
the generalized Schrodinger equation, establish a correspondence between the spaces of solutions to the initial and 
transformed equations. We show how the procedure can be used for generating families of Hamilltonians with a prede- 
termined spectrum, removing or adding new bound states. Darboux transformations for the generalized Schrodinger 
equation are considered in the differential and integral forms. It should be noted, the integral transformations can be 
important for generation of new potentials completely isospectral to a given initial one without changing the spectrum 
and for construction of Hamilltonians differing by one bound state. A correspondence between the first- and second- 
order differential Darboux transformations and the integral ones is established. To explain our findings, we consider 
a few different applications of our approach and show how it is possible to generate completely isospectral potentials 
or to construct potentials with addition or moving of bound states from the spectrum of the initial Hamiltonian. We 
analyze the influence of the distances between the energy levels on the shape of potentials, which can be very complex. 
We construct asymmetric double well and even triple well potentials for the effective mass Schrodinger equation. 

The paper is organized as follows. Section 2 is devoted to the generalized Darboux transformation of the first-order 
and corresponding supersymmetry formalism. In Section 3 we construct chains of generalized Darboux transformations 
by iteration of first-order Darboux transformations. Furthermore, we derive relations for potentials and solutions, 
obtained within the higher order Darboux transformations, in terms of potentials and solutions of the initial equations 
with no use of ones for the intermediate equations. In Section 4 we derive transformations in an integral form and 
establish relationship with differential Darboux transformations, afterwards we consider totally isospectral potentials. 
In Section 5 we illustrate our generalized transformations by concrete examples. 



II. THE INTERTWINING RELATIONS AND SUPERSYMMETRY 

As it is well known, Darboux transformations (or supersymmetry) in nonrelativistic quantum mechanics allows one 
to produce Hamiltonians whose spectra can differ in one bound state. Recently [s^, we have considered the first-order 
Darboux transformation for the Schrodinger equation with a position-dependent mass and weighted energy 
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Here m{x) stands for the particle's effective mass, q{x) and v{x) denote the potentials, (j){x) is the wave function 
and E denotes the real-valued energy and we use atomic units. In this section we will bricfiy describe the method of 
intertwining operators of the first-order, then we will analyze the supersymmetry and the nth order transformations. 



The first-order Darboux transformations 



Consider two generalized Schrodinger equations 
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where Hamiltonians Ti. and Ti. differ only in potentials v and v. 
Schrodinger equation two one-dimensional Hamiltonians % — —dx 
(j) and (j) are related by means of differential operator C 

6 ^ C6 . 



(2) 
(3) 



In the conventional intertwining technique for the 
E -f I) and % = —dxx + v and corresponding solutions 



(4) 
(5) 



As it is known, the method of intertwining differential operators provides an universal approach to generating new 
exactly solvable equations [isl [l^ . We expand the intertwining relations on our generalized equation in order to 
determine the intertwining operator £, a new potential and solutions of transformed equation ([3]) assuming that the 
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solutions to the initial equation ([2]) are known. We seek for the intertwiner C in the form of a linear, first-order 
differential operator 



C = A{x) + B{x)d/dx, 



(6) 



where the coefficients A and B are to be determined. To this end we insert ([6]) and the explicit form of the Hamiltonians 
H and T-l into the intertwining relation ^ and apply it to the solution (f> of Assuming the linear independence 
of derivative operators of different order d'^/dx'^, k = 0,1, 2, 3, we obtain the following system of equations for A, B 
and V 
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where the prime denotes differentiation with respect to x and arguments have been omitted. From ([T]) it follows 
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where /3 is an arbitrary constant of integration. Equations (jS]) and ^ allow us to determine the potential v and the 
function A. Multiplying ((51) by A and © by B and subtracting the resulting equations we obtain the equation with 
respect to A 



A" 
q m 



2A' 



B 



q m B q m B q \m 



I / l\' B' 



B 



-(- 

q \m 



A-B 



-i- 

q ym 



A' = 0. 



(11) 



In order to solve this equation, we introduce a new auxiliary function K defined by = A/B. Taking (jlOp into 
account, after simplification we arrive at the following nonlinear equation 
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from which the generalized Riccati equation follows 

— i-K + if 2 
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(12) 



where A is an arbitrary integration constant. Equation (jl2p turns into the conventional Riccati equation in particular 
cases of TO = const and q — const. With B from (|10|) and A — K/B the potential v can be expressed from ([8]) in 
terms of K and known potentials v, to and q 



V = v + 2 
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The potential v will be finally determined after finding the function K from the Riccati equation (jl2p . that can be 
linearized and integrated by introducing an auxiliary function lA — U (x) 



K = 
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(14) 



Assuming that lA is twice continuously differentiable and substituting ([M)) in ([T2l) we get the equation 

vU^qXU, (15) 
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which is identical to the initial equation ^ at £ = X. Since the solutions of ([2]) with the given potentials v, q and m 
are known at all energies £, hence wc know the solution U. Once U is given, the function K is calculated via (|14p . 
which in turn determines A by means oi A ~ BK and (fTO|) 

1 U ' 



The constant of integration /3 has been taken to one. By insertion of K in ([T^ we obtain the explicit form of the 
transformed potential 
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Finally, the intcrtwiner £ and the transformation solution </> are determined from (|6|) and (O, respectively 
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From (fTB|) - (dH]) it follows that the transformed potential w, the intertwiner £ and solutions depend not only on the 
potential v but also on the additional potentials m and q. It is easy to check that all expressions reduce to the well 
known ones for the Schrodinger equation if potential functions m and q are taken to be constant. When q = const or 
m = const we get first-order Darboux transformations for position-dependent mass Schrodinger equation 34] or for 
Schrodinger equation with weighted energy (201] as particular cases of our approach. 

Solutions at the energy of transformation. Note that relation ([T8)) connects the solutions of two equations ([2]) 
and ([3]) at an arbitrary energy except for solutions at energy of transformation, £ = X. Evidently, at £ = X the action 
of Darboux transformation P7)) on the function U and on functions (j) linearly dependent to U gives us £U = 0. 
In order to obtain a solution of the transformed equation ([3]) at energy A, we replace U hy a linearly independent 
solution U 
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The limits of integration depend on the boundary conditions. The direct substitution hi into shows that lA really 
solves the generalized Schrodinger equation. The action of C on the function U gives us a solution 77 of the transformed 
equation Q at energy A 
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By using the generalized Liouville formula ([T51) once more, one can get a second solution 77 of ([21) at energy A. 
this U is replaced by r] in and with (^0)) we get 
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Hence, the knowledge of all solutions of the initial equation © provides the knowledge of all solutions of the trans- 
formed equation ([3]) , including the solutions at energy of transformation. Notice that if U replies to the bound state 
of "H, then the function 77 defined by (|20p at the energy of transformation A cannot be normalized. Such bound state 
is excluded from the spectrum of %. Therefore, the Hamiltonians % and T-L are isospectral except for the bound state 
with energy A, which is removed from the spectrum of %. 



B. Supersymmetry 



Now we consider the generalization of the supersymmetric formalism (SUSY) on the generalized Schrodinger equa- 
tion ([1]) and show how one can construct a Hamiltonian with an additional bound state by using supersymmetry. As 
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it is well known, the SUSY algebra provides a relation between superpartner Hamiltonians, which can be presented in 
a factorized form in terms of Darboux transformation operators L and its adjoint . In our case the scalar product 
of functions is defined by not the standard way (/, x) but with the weight of q: {f,x)q — j f*[x)l{x)x{x)- Iii ttiis 
case instead of operator £>+ = (CQ)"*" = Q^C^ it is necessary to consider the operator Z?^ = q~^D^q. Therefore, 
the operator £^ adjoint to C — A + Bd/dx is determined as £^ = q~^ {A — Bd/dx — dB/dx) q. After simplification 
we arrive at 
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and the operator £^ satisfies the intertwining relation 

nC^^C^n. (23) 
The generalized Schrodinger equations ^ and ^ can then be written as one single matrix equation in the form 
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^ 0. 



(24) 



On defining Eg = diag('H,'H) and $ 



)"^, the above matrix Schrodinger equation (j24p can be written as 

\Hs - A/] $ = 0, (25) 



where / is the 2x2 unity matrix. Similar to the case of the standard Schrodinger equation, we define two supercharge 
operators Q, as follows: 
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where £ and £^ are the operators given by ([TT]) and (|22]) . respectively. One can show that the matrix Hamiltonian 
Us-, as given in (pS)) . satisfies the following conditions: 



{Q,Q} - {Q^g^} - , 
{Q,H,\ = [i7„gt] =0, 



(27) 
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where {•, •} and [•, •] are the anticommutator and the commutator, respectively. The relations ([27l) are trivially fulfilled, 
because the matrices in (|26p are nilpotent. It is easily seen that equations (j28p are equivalent to the intertwining ones 
O and 

Now, let us consider the complementing relations of the supersymmetric algebra, that is, the anticommutators 
{g, g^} and {g^ g}. For this, we calculate the operators C) C and ££^, and consider the connections of them with 

our Hamiltonians and l-l. By using ([T7|) and (|22|) . we arrive after some algebraic transformations at 
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We express the potential v from the Riccati equation ((12]) in the form 

v = — i-K' + K'^)-(—\ K + qX . 
m \m J 

Using this in (|13|) , we get the following representation for transformed potential 
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One can easily see that the potential difference is determined as 

If d K 
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On further employing ([3T|) and (|32|) . the formulae ([29)) and ([30]) can be rewritten as 



£^£ - 

As can be easily seen, the anticommutation relation is 
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Note that as soon as the initial Hamiltonian H is presented in the factorized form ([37|) , one can get its supersymmetric 
partner in a factorized form ([38]), too. Indeed, multiplying equation (|37p from the left by C and taking into account 
the intertwining relation ()23p we get 



L {C^C + X)<t> = {C + X)C <t> = 'H C 



(39) 



It means, that n = CC'^ + X. 

In summary, we obtained the explicit forms of the supersymmetric partner Hamiltonians % and %. The Hamil- 
tonians ([37| and ([38| are compatible with their definitions ([2]) and ([3]), respectively, if the transformed potentials v 
and V are given by ([3T|) and ([5^ . Finally, taking the difference of the factorized Hamiltonians ([57| and (155)) gives 
the potential difference (|13p that we obtained for our Darboux transformation. Hence, the Darboux transformation 
is equivalent to the supersymmetry formalism. 

The intertwining relation (|23p means that the operator £t is also the transformation operator and realizes the 
transformation from the solutions of ^ to the solutions of Evidently, one can interchange the role of the initial 
generalized Schrodinger equation ([2]) and its transformed counterpart ([3]). To this end, let us express the operators C 
and in terms of functions rj given in (1201) . which are solutions to the equation ([3]) at the energy of transformation 
A. First, we rewrite K by using U from the relation ([20)) 
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Using this in ([TTI and ([22)) . we obtain after simplifications 
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Obviously, the function 77 is also a transformation function. Notice, L^rj — 0, meaning that rj belongs to the kernel of the 
operator £^ . As one can see from ([101) and ([^T)) , the application of the operator £t to the second linearly independent 
solution fj of equation ([3]) gives back the solutions U of the initial problem at the energy of transformation. Indeed, 
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Hence, the solution at the energy of transformation A takes the form 
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and, in principle, can reply to the new bound state. The second linearly independent solution of ([2|) at energy A 
can be written in terms of 77 as follows: 



\ q r] 
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Introducing the function K = K{x) hy K = rj' /-q and taking into account -^^=f {\\ — 57^ 
expressions (|40p for operators and C can be rewritten as 



(42) 



^ /X the 

q ax V m ' 



C 



1 



K 



c 



1 



d 



/q m \ dx J ^q m \dx 

Using (|33p the potential v can be expressed in terms of v as follows: 

d ( k \ d /I d ^ 



K 



V = V ■ 



dx 



and corresponding solution (j) are given by 



dx \q dx 



d 

dx 



idfV^ 

q dx \ m 



K 



(43) 



(44) 



(45) 



Thus, the function rj becomes a transformation function for the operator C\ which performs the transformation from 
the potential v to the potential v and from the solutions of ^ to the solutions of ([2]). If within the first procedure 
(|17p - (|18p we constructed the potential v with one bound state removed, now we can construct the potential v with an 
additional bound state. Note that we have established a one-to-one correspondence between the spaces of solutions of 
equations ([2]) and ([3|). The operators C and realize this correspondence for any £ ^ \. If f = A, the correspondence 
is ensured by mapping r/ < — > 11 and r\ i — > hi. 

In particular cases our generalized Darboux transformations are reduced correctly to the known expressions. In 
the case with a constant weighted energy potential, e. g. q(x) = 1, from our supersymmetry approach we get the 
supersymmetry for effective mass Schrodinger equation [29l, |34!|. In the case with constant mass m{x) — mo from our 
approach we obtain the supersymmetry for Schrodinger equation with weighted energy 20]. Finally, if m(x) = mg and 
q{x) = 1, our expressions of supersymmetric algebra are correctly reduced to the conventional ones for the standard 
Schrodinger equation (see, e.g. [Sl). 



III. HIGHER ORDER DARBOUX TRANSFORMATIONS 



In this section by considering iterative applications of first-order Darboux transformations n times we obtain the 
nth order Darboux transformation and show that the nth order Darboux transformation can be expressed in terms 
of solutions of an initial equation, with no use of the solutions to intermediate equations. 



A. Chain of Darboux transformations 



Now we consider iteration of one-step transformations we have obtained in the previous section, in order to construct 
higher order Darboux transformations. To this end we show that the intertwining operator £ of the nth order can be 
obtained from a sequence of n first-order Darboux transformations, like conventional Schrodinger equation, and create 
a chain of exactly solvable Hamiltonians ^1,^2, ■■■^'Hn- Let us define the second-order Darboux transformation as a 
sequence of two first-order Darboux transformations Ci and C2 

C = C2C1 , (46) 

where £1 is actually L given in ([T7|) . Let Ui be an auxiliary solution of ^ at energy Ai. Then we have 



8 



The operator £2 is determined as follows: 
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The function xi is obtained by means of the first-order Darboux transformation ()47p . applied to an auxiliary solution 
U2 of equation ^ at energy A2 
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Clearly, xi is the solution of the transformed equation with the potential vi = v, defined as in (jl3p . and xi can 
be taken as a new transformation function for the Hamiltonian Hi — H to generate a new potential 
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The function ipi — (j), defined as in p^ . is the solution of equation ([3]) with the Hamiltonian Hi = H. In summary, 
the action of the 2nd-order operator on solutions (p of the generalized equation ^ leads to solutions 02 of 
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given by 



— C (f) = C2 Ci 
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Iterating the procedure n times in regard to the given operator % leads to the operator H„, which satisfies the 
intertwining relation 



The transformed potentials Vn satisfy the following recursion relation 
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where C is the nth order operator: 
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Thus, the chain of n first-order Darboux transformations results in a chain of exactly solvable Hamiltonians % — > 
T-Li — > ... — > T-Ln- When C is the nth order differential operator and the intertwining relation ((54| is valid, the so-called 
nth order supersymmetry arises, like for the ordinary Schrodinger equation [To| . 



B. Darboux transformation of the nth order 



Now we show that iteration procedure of the Ist-order Darboux transformations in (|55p - (j57p can be removed and 
transformed potentials u„ and solutions 0„ can be expressed in terms of the initial potentials v, m and q and the 
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family {Uj}, j = 1, 2, n of auxiliary solutions of the initial equation ([2]) at energies \j. Consider now the 2nd order 
transformation in detail. Using the explicit expression for vi which appears in the first-order Darboux transformation 
([T3l) . we present the formula (|50|) for the potential V2 as 



V2=V + 2 
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where Wi,2 — U1U2 —U1U2 is the Wronskian of the functions Ui and U2- Substituting (l59l) into the formula (|48)) for 

K2, we get 
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With the last expression after some manipulations, the new potential V2 can be expressed as 
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and finally in the form 
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Find now the corresponding functions (j)2. To this end let us transform the relation (|5ip. By analogy with xi the 
function (pi can be written in terms of the Wronskian Wi,£ — Ui(p' — Uic/)' 
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Let us now calculate the derivative of 1 
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Making use of the last expression and the relation ([60]) for K2 in ([5T|) . we obtain, after some simplification, the 
solutions as follows 
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It is easily seen from (j63p and (j65p that due to the 2nd order Darboux transformations, the potential and solutions 
are completely expressed in terms of the known potential functions, w, m and q and the solutions Ui,U2,(f>{£) of the 
initial equation, with no use of the solutions to the intermediate equation with the potential vi (x) . 

Clearly, for the next transformation step to be made, one should take a new transformation function X2, that 
corresponds to the potential V2 at energy A3. The solution X2 can be obtained by applying the operator C — C2C1 to 
the transformation solution at energy A3, that is 



X2 = C2C1U3. 
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According to (1651) . the solution xi can be written as 



X2 



1 



23 



q m Wi2 

and can be used to produce a new transformed operator £3, given by 



(66) 



C3 



1 



/q m \dx 



d 



A3 = 

X2 



for generating a new potential V3 with corresponding solutions (/)3 and so on according to ([55)1 - ([571) . In this way we 
can express the transformed potentials Vn of any order in terms of the initial potentials v, q, the effective mass m and 
the family of auxiliary solutions Uj, j = 1,2,, ...n of the initial equation ^ at energies Xj, which different from each 
other: 



Vn = V + 2 



q d K 



m dx Jq m 



q d 
n \l—-r 

m ax 



1 d 

q dx 



K = Ki+K2 + 



(67) 



In the case of n-order transformation over the initial potential v equation (|67p gives us the n-SUSY partner potential 
Vn- This construction enables us to generate a family of new Hamiltonians of any order and corresponding solutions 
directly from the initial Hamiltonian and solutions without generating intermediate Hamiltonians. In this section we 
constructed nth order Darboux transformations by using a chain of first-order Darboux transformations and shown 
that the nth order Darboux transformations are equivalent to the resulting action of the iterative method. 

IV. THE INTEGRAL FORM OF DARBOUX TRANSFORMATIONS. 

In this section the generalized Darboux transformations are represented in the integral form and applied to con- 
struction of Hamiltonians with the same spectrum as the initial one (totally isospectral Hamiltonians) , and differing 
by one bound state and by two bound states. 



A. First- and second-order integral transformations 



At the beginning consider the first-order transformation. Multiplying both sides of equation ([2]) for (j) by Ui and 
subtracting from the obtained expression the equation ([S]) for Ui multiplied by (j), we arrive at 



d /Wl^£ 



dx \ m 

The last expression can be easily integrated: 



(Ai-£) qUi 



Wi,£ = m 



{\i-E)j q{x')Ui{x')(f>{x')dx' + , 



(68) 



(69) 



where C is a constant of integration. Inserting the integration result into the formula (j64p for 01, we arrive at the 
integral form for the 1st order transformed solutions: 



C+{Xi~£) / q{x')Ui{x')(f,{x')dx' 



m 1 
qlTi 

The auxiliary solutions xi, taking at the energy £ — X2, can be written as 



XI - (c\ + (Ai - A2) q{x')U,{x')U2{x')dx'^ 



(70) 



(71) 



where it was used Wi_2 — ni ((Ai — A2) q 14x142 dx' + Ci). The integration limits depend on the boundary condi- 
tions. 

Now consider tlie second-order transformation. By analogy, using 

W^xi,0i="^ f(A2-f) / q{x')xi{x')Mx')dx' + C 
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in ([5T|) . we shall get the integral form for the 2nd order transformed solutions presented in terms of the 1st order 
solutions <pi and Xi 



»2 = 



/q m xi 



y^-^ {c + (A2 - E) qix')xi{x')M^')dx''^ . (72) 

Evidently, the solutions 02 can be expressed directly in terms of the solutions to the initial equation For 
this transform the expression (I55|) for 02 having regard to equation ^ for (f> and Ui, using ([M]) and ;^[M^i,2]' — 

[m]'^i-2 obtain 



Wl,2 ' 



m[Ci + (Ai - A2) / q{x')Ui{x')U2{x')dx'] 
^ , , ^ w ^2 [C + (Ai - g) (xO0(xOdx'] 

c,+rq{x')U,{x')U2{x')dx' ' ^^"^^ 

where we have introduced ci = Ci/(Ai — A2). Using the integral presentation of Wronskian in (|63|) . the transformed 
potential W2 can be written as 

^ ^ _ 2 /TA (J: g^2^^i A .74^ 

" \l mdx ci+ f dx' q{x')Ui{x')U2{x')) 

Thus, we get the integral form of the first- and second-order Darboux transformations for the potentials and solutions. 
Note, the spectrum of Hamiltonian 7^2 with the potential V2 differs from the spectrum of the initial Hamiltonian H 
by two bound states Ai and A2. 

The formulae (j74p and (|73p can be rewritten in a more simple form 



V2ix)^v{x)~2J^J-l \ ^ Kix,x)] , (75) 
y m,{x) dx \^q{x)m{x) J 



K{x,x')<j){x')dx' (76) 



with the operator kernel K{x, x') determined as 



Kix,x')^ ^2Wg(^0^i(^0 (77) 

^ ' ci + J dx' q{x')Ui{x')U2{x') 

In (|76l) the constant C, connected with the Wronskian Wi^g, is chosen to be zero. It is interesting to note, that at 
q{x) = const and m(x) = const the formulae (j75p and (j76p look like the inverse problem ones for potentials and 
solutions [Ui, 17] obtained with the degenerate kernels K{x,x') except for the form of transformation operator (|77p . 
The operator kernel (|77p differs from the inverse problem kernel K{x, x') not only by q{x) ^ const but by auxiliary 
functions Ui{x) and U2{x), which correspond different energies Ai ^ A2. In the next section we shall consider the case 
when the 2nd order Darboux transformations allow to change the spectrum of a given Hamiltonian on one bound 
state that corresponds to Ui{x) = U2{x). 



B. Hamiltonians 14.2 differing by one bound state and completely isospectral Hamiltonians 

The first-order supersymmetry, above considered in the differential and integral forms, gives us opportunities 
to construct isospectral Hamiltonians differing by one bound state. It is interesting to note, if we use Darboux 
transformation in its integral form, then we directly from (1701) obtain the solution of the partner equation (jSj at 
energy of transformation £ = \i which with an accuracy of an arbitrary constant coincides with (j20l) 




Now we show how using double Darboux transformations to generate Hamiltonians 7^2 1 the spectrum of which 
differ from the spectrum of the initial Hamiltonian % by one bound state and how to construct a family of isospectral 
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Hamiltonians 'H2^ the spectrum of which completely coincide with the spectrum of %. For the 2-nd order transfor- 



U2 



mation we used the transformation function Xi{^) obtained within the first step, xi — ^il^2 — ~ TJ^ 

with Ai 7^ A2. But for A2 = Ai we have hl2 = lAi and it leads to xi = i^il^2 = 0. For the second transformation one 
can use the function xi constructed by means of a linear combination of the solutions r] and r): xi — ^iV + ^2^- For 
our aim we take a linear combination as follows 



Xi = C177 + 77 = \ (^ci + j dx' q{x')Ul{x')^ , 



(78) 



By analogy with (|FT|) we calculate K — Ki + K2 with Ki = —Ui/Ui and K2 ~ —Xi Ixi- After simplification we 
have 



_d_ 

dx 



(ci + J"' dx'q 



and 



Kix)^--[lnJ-, 



[ci + f dx'q Ul) 



(ci + fdx'qUf) ■ 



(79) 



Plugging the last expression into the formula (|58p which defines the potential, after some transformations we arrive 
at 



V2 



q d 



mdx \\ m ci+ J"" dx' q{x')Ul{x') , 
With xi defined by ([75]) and represented by its integral form ([70)1 the relation ((7^ leads to 



(80) 



^2 ^ 0(A. - 8) - ^^^^.J^^2^,^^, {C + (A, - 8) f qix')U.ix')^ix')dx) . (81) 



It should be noted, that in difference from the differential approach, the relations ((80)) and ([8T|) for the new potential 
V2 and the solution 02 can be obtained directly from (j74p and (I73p , which replies the transformations with two bound 
states, if one takes Ui =U2- 

It is worth noting that the auxiliary function xi can be determined as follows 



Xi = ?7 + r77 



m 1 
Ilk 



1 



dx'q{x')Ul(x" 



Then the potential V2 and solutions (j)2 are rewritten in the form 



q d 



V2 ^v-2J J — 



ml 



mdx \\j m 1 + r/"' dx' q{x')Ul{x') 



(82) 



(83) 



{\i-£)4>- 



TUi 

1 + r f q{x')Ul{x')dx' 



C+{\i-£) j q{x')Ui{x')(j){x')dx' 



(84) 



In principle, the formulae (|82l) - (|84|) coincide with (f78| . (|80|) . (|8T|) at F = 1/ci, but they are more suitable for physical 
applications. Here, the constant F plays a role of a normalization constant of the new bound state Ai provided the 
other spectral characteristics coincide. 

By analogy with the differential Darboux transformations, the solution of the generalized equation ([T]) with the 
potential (j83l) at energy of transformation £ = Xi can be achieved by means of operator C2 acting on the solution rj 
from (I20p . obtained within the first transformation step 



m = C2V = 



1 



d 

dx 



Xi 



m 1 
7^ 



(85) 
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After transformations (1851) with an account of ((82|) we get 

" ~l + Tj^dx'qix')UUx')' ^^^^ 

Note, the solution 772 can be directly obtained from (|M| at £ = Ai. Obtaining solutions at energies of transformation 
from the general formulae at arbitrary energies is one of advantages of integral transformations. The relations for 
potential (I55|) and solutions (IM|) can be rewritten in the form ([75]) and (|7^ with the operator kernel K{x,x') 
determined as 

K{x,x') = 



TUi{x)q{x') Ui{x') 
1 + r/'' dx' q{x')Ul{x') 



At m{x) = const = mo and q{x) = const our generalized expressions are correctly reduced to the integral equations 
of inverse problem for the standard Schrodinger equation with the degenerate kernel of transformation (see e.g. 

SElii)- 

It is worth mentioning that the double Darboux transformation with eliminating or adding one bound state at 
an arbitrary energy allows one to avoid the problems with singularities of the transformation kernel K{x) and, as 
consequence, to avoid the problems with singularities of the constructed potentials and solutions. We assume, that 
771 (a;) and q{x) do not lead to the additional singularities on the tested interval. Let us compare the formulae for the 
first-order transformation (fT7|) . (IT5|) and for the second-order one (1751) . (ESI)- One can see, if the first step procedure 
is based on any arbitrary solution U of the generalized equation ([2]) , rather than the ground state wave function, the 
superpotential K = —U'/U becomes singular. Singularities are localized at the zeros of exited wave functions. It 
leads to singularities of constructed potentials and solutions. Unlike the 1st order transformations the transformation 
kernel K{x) of the 2nd order (ffO]) has no singularities at zeros of exited wave functions. It means that one can 
make transformations on an arbitrary bound state ( not only on the ground state) and construct the potentials and 
corresponding solutions without additional singularities. 

It is interesting to note that by using the double Darboux transformations (j83p - (|86p we can construct new potentials 
V2 without changing the spectrum of the initial potential v, i.e. iully isospectral potentials. Indeed, if the bound state 
Ai belongs to the spectrum of the initial Hamiltonian H and T = N2 — is a difference between the normalization 
constants of the bound state Ai for H2 and H, the formulae (|83p - ([86]) give us a family of isospectral potentials and 
corresponding solutions, since the normalization constants can be chosen arbitrary. In quantum mechanics potentials 
whose spectra coincide and differ only in the normalizations factors A'2 and N of bound states are called phase- 
equivalent potentials. Note, the phase-equivalent potentials have a different shape. They can be more deeper and 
narrow or more shallow and wider and possess the same spectral data, except for normalization constants. 

If we assume, the transformation function lAi{x) to be taken at the energy of the bound state, which we would like 
to add to the initial spectrum, and F = 7V| is the corresponding normalization constant, then the formulae (I83p give 
us the possibility to construct a family of two-parametric potentials with a new bound state Ai and an arbitrary F, 
whereas the other spectral characteristics of the spectra produced by the potentials ^2(2:^) and v{x), coincide. Note, 
constructed potentials V2{x) from this family are isospectral among themselves, since posses the coinciding spectra 
and differ only by the normalization constants. 



V. APPLICATION 



The 1-st example. As an illustrative example we present the transformed potential and solutions corresponding 
to the first, second and third-order Darboux transformations. We start with the generalized Schrodinger equation ([T]) 
with the repulsive coulomb potential v{x) = 1/(4 x) 



1 



{x) + ^<j){x) = q{x)£(j){x) 
Ax 



dx \m{x) ) dx 

where we choose the effective mass as m — 1/x and q = x. The general solution of this equation can be written as 



(87) 



</)(x) 



Ci sin(fca;) ^ C2 cos(fca;) 



ks/x 



ky/x 



(88) 



Now we would like to generate potentials with one bound state at the energy £1 = —k\ and obtain corresponding 
solutions by a first-order supersymmetry transformation ()44|) and (I45p applied to a special case 77 = 771 of the general 
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solution (ESI) 



Vi = 



C cosh(Kia;) 



We obtain the transformation operator Ki — rj'i/rji in the form 



Ki = h Ki tanhKi x, 

2x 



the potential vi and corresponding solutions (f>i at £ ^ £i 



viix) — 2xKi 1 — tanh (kix)) = ^ . 

^ ' ix ^ ^ ^ ' ix cosh2(Ki x)' 



(89) 







d 1 








dx 




dx 2x 



Ki tanh Ki X 

The solution at the energy of transformation £ = Xi = — Kj is defined in accordance with (|4T|) as 

y q Vi vxcosh^Ki X) 



(90) 



(91) 



we obtain 



and corresponds to the bound state. Note, by varying (j) in (|90|) we recover all solutions of ([TJ with the transformed 
potential (j89p and with m = 1/x and q — x. In the case when is chosen to be (t>{x) — fc^"^^ 

C cos(/cx) C Ki tanh(Kia;) sin(fca;) 



y/x k^/x 

If (f) is chosen as </) = '^^k!^^^^^ ^'^^ following partial solutions 

,+ , ., , , C exp{±ikx) 
(pi = (^ZK + Ki tanh(Kix)) — j= 

/Cy X 



(92) 



(93) 



Thus, we presented the simplest example of exactly solvable problem for the generalized equation ([T]) with q{x) — x, 
m{x) = 1/x and with a real potential (|89p . which is singular at zero. The potentials, obtained at different energies of 
transformation, are depicted in Fig. la. 

Now employing Darboux transformations of the second-order, we shall construct potentials and solutions of the 
generalized equation with two bound states. We define the auxiliary functions rji and 772 as follows 

cosh(Kix) sinh(K2a;) 

m = — , , m - 

For the second step we have K = K\ + Ki , where 



Ai = — , if 2 = — , Xi = -??2 + ^I'ya , K = -—\\vl 

771 xi m \ / ax \ y/q m 

The transformed potential V2 having two bound states at energies Ai = — and A2 — —1^2 can be written as 

2^9 d ( K 



V2 =V - 



m dx \ Jq ml y m ax 



ld_ 
q dx 



Finally, for our choice of v, m and q we obtain 



(94) 



(95) 
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FIG. 1: (a) Potentials 111(2;) corresponding to one bound state created at different energies. (b) Completely isospectral 
potentials with bound state energy £1 = —16.0, dashed line corresponds to the initial potential vi, solid lines correspond to its 
isospectral potentials, (c) Potentials «„, n = 1,2,3 having one, two and three bound states, respectively. 



where Wi^2 — W{rii, 772) — (^'^2 cosh(Kia;) cosh(K2a:) — ki sinh(K2a;) sinh(Kix) j and the corresponding solutions 

are 



cosh(Kix) 1 dx 



d d{\nWi,2 



dx 



k^E ■ We put 



where Wi^s — ^ kj/Kx cosh(Kia;) cos(fcx) — ki sinh(Kia;) sin(fcx)^ if <f>{x) is chosen as (l){x 

Ki < K2- The potential having two bound states are depicted in Fig.lc. 

By using ([SS]) one can construct the potential V3 for the generalized Schrodinger equation ([3]) having three bound 

states 



V3 



13 

Ax 



^2 

2a;^ In 2,3 



(96) 



where the auxiliary functions rji, 772 and 773 determined as 



m 



cosh(Kia;) 
^/kTx 



sinh(K2x) 
yJWx 



m 



cosh(K3x) 

y/K3 X 



give us the Wronskian Wi^2,3 

1 



2,3 



cosh(Kia;) cosh(K2a;)K2 sinh(K3a;)K3 — cosh(Kix) sinh(K2a;)K2 cosh(K3a::)K3 

— sinh(Kia:)Ki sinh(K2x) sinh(K3a;)K3 + sinh(Kia;)Ki sinh(K3x) sinh(K2a;)K2 + 
+ cosh(Kia;) sinh(K2a:)Kj cosh(K3x)K3 — cosh(Kia;)Kj cosh(K2x)K2 sinh(K3a::) . 



(97) 



As an illustrative example we present the potentials f 1, obtained at the energy of transformation £1 — —4, V2 obtained 
at the energies £1 = —4, £2 ^ —16 and U3 calculated at the energies of transformation £1 = —4, £2 = —16, £3 = —25. 
They are depicted in Fig.lc. 

Now we are going to construct isospectral potentials. As an initial potential we take the potential vi, obtained 
within the first-order intertwining (|89)) . Using (l83l) and changing v by vi from (|89| and Ui by U from ((9T|) after 
simplification we obtain a family of potentials with the same eigen-value and different values T 



V3{x) 



1 2xKi d , , 

5 2x — in [1 + 1 tanh ki x\ 

Ax cosh (ki x) dx 



(98) 



The initial potential v\ and its isospectral potentials (|98l) are presented in Fig. lb. 

The 2-nd example. As the following example let us consider creation of new bound states for effective mass 
Schrodinger equation 



1 



dx \m{x) J dx 



(j){x) + v(x)4>{x) — £ 4>(x) 



(99) 
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that corresponds to the generalized equation ([T]) with q{x) = 1. We choose the effective mass in the form m{x) = o? jx^ ^ 
the initial potential v(x) — and start with the equation 



d 

dx 



1 



m(x) 



d_ 

dx 



(j){x) — E (j>{x). 



(100) 



The general solution of eauation (jlOOp can be written as 



— — = [ci sin(az^ ln(x)) + C2 cos(q:z^ ln(a:))] , 



(101) 



where ci, C2 are free constants and i^^ — (— 1 + 4Q!^fc^)/4a^. Recently in [sj] we have constructed the potentials for 
effective mass Schrodinger equation ([M)) with creation of one and two bound states without investigation of potential 
forms. Here we would like to apply our technique to construction of double-well and even triple-well potentials with 
creation of two and three bound states, to generation of completely isospectral potentials and to investigation of the 
influence of position dependent mass on the form of constructed potentials. 

Construction of completely isospectral potentials. As an initial potential we take the potential vi , obtained in (s^ 
within the 1-st order intertwining 



cosh (a7ln(x)) 



(102) 



It can be easily obtained from (|44)) at q(x)=l with a particular solution of 771 given as r]i{x) — cosh (071 ln(a;)) 
The solution at energy of transformation £1 — —k^ is 



Vm*{x) 

mix) 



fx cosh (a7 ln(x) ) 



(103) 



and corresponds to the bound state. Using (1551) and replacing v by v\ from (|102p and U\ by U from (|103p after 
simplification we obtain two-parametric family of isospectral potentials 



27' 



cosh (a7ln(.T)) 



where 



2x d 

a dx 



X d 



' dx 



InP 



(104) 



p = 1 + r 



cosh (a7ln(x)) 



-dx . 



All these potentials V3 posses a single bound state each with the same energy £1 = — as the initial potential vi, 
as well as the normalization constants including in T can be chosen arbitrary. In Fig. 2a we have plotted the initial 
potential calculated by the formula (jl02p and its strictly isospectral potentials, calculated by (|104[) at different F . 
Influence of distance between levels and effective mass on the form of potentials. By using second-order intertwining 
let us construct a potential with creation of two bound states. For this we define the auxiliary transformation functions 
as follows 



cosh (a7i ln(x) ) , 772 [x) 



sinh (q!72 ln(x)) , 



(105) 



where 



(1 + 4q! Kf)/4Q! , i = 1,2. We put ki < K2. The potential V2, obtained within the second-order 



Darboux transformation, can be expressed from (1941) with q{x) = 1 



V2 



2 d 
m dx 



1 d W: 



12 



VF12 dx 



2 d^ 



m dx"^ 



(106) 



This formula coincides with expression obtained in [3J] for effective mass Schrodinger equation. For our choice of 
m,{x) the last term vanishes and the potential is written as 



V2^ 



2x d 
a dx 



- ^\nW^2 
a dx 



(107) 
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FIG. 2: (a) Completely isospectral potentials with bound state energy £i = —2.0, dashed line cor- 

responds to the initial potential vi, solid lines correspond to isospectral partner potentials. (b) The 
change in form of potentials V2{x) with two bound states as the levels approach each other: 
1. — >■ = —2.0, £2 = —6.0, 2. £1 = —2.0, £2 — —3.75. (c) The influence of m{x) on the behavior of transformed 
potentials V2{x) with two bound states £1 = —2.0, £2 ~ —3.75. 

2.5 5.0 7.5 10. □ 12.5 l^i 2.5 5.0 7.5 10.0 12.5 2.5 5.0 7.5 10.0 12.5 
□ K J □ I '! O k x 



(a) 




FIG. 3: The change in form of potentials Vs{x) with tree bound states as the levels come close to each other: 
(a) 1. ^ £i = -1.0, £2 = -3.75, £3 = -6.25, 2. ^ £1 ^ -1.0, £2 = -3.75, £3 = -5.0; (6) £1 = -2.0, 
£2 = -4.75, £3 = -6.0; (c) £1 = -2.0, £2 = -3.75, £3 = -5.0. 



where Wronskian is determined as 

2 

W12 = — ^72 cosh(a72 ln(x)) cosh(a7i ln(x)) — 71 sinh(a7i Inx) sinh(a72 ln(a;))^ . 

The potentials having two bound states are presented in Fig. 2 b,c. The graphs in (Fig. 2 b) depict the forms of 
constructed potentials in dependence from the distance between energy levels. One can see if the levels are close to 
each other, we construct asymmetric double well potentials (curves in Fig. 2 2b) another, we construct asymmetric 
potentials (Fig. 2 lb). Note, double well potentials have attracted some attention over the last years (see, e.g., [iol - li^ ). 
Asymmetric double well potentials for the ordinary Schrodinger equation were investigated in (40| with introducing 
a special parameter of asymmetry. In our case asymmetry in forms is a consequence of the position-dependent mass 
m{x) that is singular at zero. The influence of rn(a;) on the form of constructed potential v{x) is demonstrated in 
Fig. 2 c,b. Our analysis shows that the larger a, the shallow and narrow constructed potential (see Fig. 2 b, curve 2 
and Fig. 2 c). Another words, increasing m{x) leads to decreasing potential v{x). 

The next considered example illustrates the possibility to construct potentials having three bound states, for all 
that we can generate double well potentials and even triple well potentials. Employing the third-order Darboux 
transformations (j67p with q{x) — 1, the potential V3 can be written as 



V3 



2x 


d 


X d 


a 


dx 


a dx 



■InWi 
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(108) 



The Wronskian Wi,2,3 is determined by the auxiliary functions r]i{x), 772(2:), Vai^)^ where ?7i(x), 772(2;), is defined in 
(|105l) and 773 is given as 773 (x) = cosh (073 ln(a;)) . The potentials V3 calculated by the formula (|108l) are plotted 
in Fig. 3. As in the previous case with two bound states, the forms of constructed potentials depend on the space 
between energy levels. We can see if levels are sufficiently distant from one another, we construct simple asymmetric 
potentials presented in Fig. 3 a (curve 1), if two levels out of three are close to each other, we construct asymmetric 
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double well potentials (Fig. 3 a, curve 2), if three levels are close to each other, we construct asymmetric triple well 
potentials (Fig. 3 b, and Fig. 3 c). As a final remark, let us note that different distances between levels give us different 
shapes of potentials. It can be very important for construction and investigation of quantum systems with needed 
spectral properties, e.g. in nanoelectronics (i^ . 



Conclusion 



By application of the intertwining operator technique to generalized Schrodinger equation with position dependent 
mass and with energy dependent potentials, Darboux transformations of an arbitrary order have been constructed. It 
has been shown that nth order Darboux transformation is equivalent to the resulting action of a chain of first-order 
Darboux transformations. Our generalized Darboux transformations comprises the position-dependent effective mass 
case and the case of linearly energy- dependent potentials, as well as the conventional case of Schrodinger equation. 
The integral Darboux transformation method has been elaborated for the generalized Schrodinger equation. An inter- 
relation has been found between the differential and integral transformations. The integral Darboux transformations 
have been applied to generation of isospectral Hamiltonians differing by one and by two bound states from the spec- 
trum of the initial one. It has been shown how to produce completely isospectral Hamiltonians to a given initial one. 
On concrete examples it has been demonstrated how to apply the Darboux transformation technique for modeling 
quantum well potentials with the given spectrum. Hamiltonians with different number of levels have been produced 
and the influence of the distance between levels on the shape of constructed potentials has been investigated, in 
particular, asymmetric double well and triple well potentials have been built. The influence of the position-dependent 
mass on the behaviour of constructed potentials has been studied, too. 
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